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According to the second law of thermodynamics the total entropy of a sys- 
tem is increased during almost any dynamical process. The positivity of the 
specific heat implies that the entropy increase is associated with heating. This 
is generally true both at the single particle level, like in the Fermi acceleration 
mechanism of charged particles reflected by magnetic mirrors, and for complex 
systems in everyday devices. Notable exceptions are known in noninteracting 
systems of particles moving in periodic potentials. Here the phenomenon of 
dynamical localization can prevent heating beyond certain threshold. The dy- 
namical localization is known to occur both at classical (Fermi-Ulam model) 
and at quantum levels (kicked rotor). However, it was believed that driven er- 
godic systems will always heat without bound. Here, on the contrary, we report 
strong evidence of dynamical localization transition in periodically driven er- 
godic systems in the thermodynamic limit. This phenomenon is reminiscent of 
many-body localization in energy space. 

PACS numbers: 05.30.rt: Quantum phase transitions, 05.50.+q: Lattice theory and statistics, 
05.45.-a: Nonlinear dynamics and chaos 

I. INTRODUCTION 

The second law of thermodynamics puts very strong constraints on dynamical processes. In 
particular, it implies that the entropy of driven closed macroscopic systems increases in time for 
almost all dynamical processes. In most cases the entropy increase is associated with heating, i.e. 
the temperature and the energy of these systems also increase. The second law of thermodynamics 
can be rigorously proven microscopically if the initial state of the system is stationary, i.e. the 
initial density matrix is diagonal in the energy basis, and the probability distribution of occupying 
different energy states is a decreasing function of energy (so called passive density matrices )i. Then, 
as a result of any dynamical process, the energy of the system can only increase or stay the same. 
This is exactly the Thompson's formulation of the second law of thermodynamics^. Similarly, with 
a single assumption of the initial state being stationary, one can prove that the properly defined 
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entropy increases or stays constant*. Even more generally without any assumption on the initial 
state, it can be proven that the entropy and energy (for positive temperature states) increase in 
time for systems subject to random forces. These systems include the famous Fermi acceleration 
problem, where a charged particle increases its energy by being repeatedly reflected from a magnetic 
mirror—, the Boltzmann's dynamics of particles experiencing repeated and uncorrelated collisions^, 
particles moving in a chaotic time-dependent cavity^ and many others. In all previous examples it 
was assumed that there were no correlations between subsequent collisions. In quantum language 
this lack of correlations is equivalent to relaxation to the diagonal ensemble^ and it is sometimes 
referred to as "energy measurements"—. 

The applicability of the second law becomes much less clear if there are correlations between 
subsequent dynamical processes, like in the case of systems which are driven periodically in time. 
In these systems one can observe the phenomenon of dynamical localization, where the energy of 
the system never exceeds a maximum bound. Examples of dynamical localization include Fermi- 
Ulam model of a particle bouncing off a periodically moving wali^, classical and quantum kicked 
roto r 11 ' 12 , and the Kapitza pendulu m 13 ' 14 . It is generally expected that this localization phe- 
nomenon is peculiar to small integrable systems. In chaotic ergodic systems the periodicity of the 
driving should not matter because such systems effectively serve as their own heat bath, which in 
turn can be viewed as a source of a random Langevin type noise. Thus, it is generally expected that 
the energy of a periodically driven ergodic system will steadily increase in time. This belief though 
is entirely based on intuition and there is no guarantee that it is generally correct. Moreover, 
using the example of the many-body localization transition, it has been convincingly argued that 
disordered interacting systems can behave non-ergodicall y 15116 and the standard assumption that 
the interacting non-integrable system can serve as its own bath can fail. Furthermore, for driven 
interacting systems it was recently noticed numerically that periodic modulation, even with a long 
period, leads to suppression of energy and entropy growth 3 - 1 ^ compared to a a similar modulation 
with a random period. However, it was not clear whether this suppression simply leads to a slower 
heating rate or to the localization. 

The purpose of this work is to argue that there is a new type of localization transition in 
periodically driven systems as a function of the driving period. We demonstrate this transition both 
analytically and numerically and argue that it is related to a breakdown of short time (Magnus) 
expansion of the evolution operator. For short periods the expansion is convergent leading to the 
effective time-independent many-body (Floquet) Hamiltonian and the energy is localized. For long 
periods this expansion breaks in the thermodynamic limit leading to the derealization transition 
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FIG. 1: Schematic representation of the Kapitza pendulum, i.e. a rigid pendulum with vertically oscillating 
point of suspension, and its phase portraits, (a) The Kapitza pendulum, (b) Undriven regime: the pendulum 
performs small oscillations around the stable lower equilibrium which are represented by the red line in the 
phase portrait, (c) Dynamical stabilization regime: the pendulum performs small oscillations around the 
stable upper equilibrium which are represented by the red line in the phase portrait. In the phase portraits 
the green lines correspond to rotations, the black lines to oscillations, the blue lines are the separatrices and 
the points represent the region of chaotic motion. 

and heating of the system to the infinite temperature. We show a strong evidence for this transition 
using examples of classical and quantum spin chains. 



Before addressing the many-particle problem we will discuss a much simpler example of a 
periodically driven system, so called Kapitza pendulum^ and show how the Magnus expansion 
can be used to derive the effective potential. The Kapitza pendulum is a classical rigid pendulum 
with a vertically oscillating point of suspension (see Fig. [I]). 

The equations of motion of the Kapitza pendulum read 



where 6 is the angle measured from the downward position (see Fig.[T]), loq = w | is the frequency of 
small oscillations and a, 7 are the amplitude and frequency of the driving of the point of suspension: 
y c = —a cos (7^). This dynamical system has an extremely rich behavior containing both regions of 
chaotic and regular motion (see Ref. and references therein). For our purposes we consider the 



II. MAGNUS EXPANSION AND KAPITZA PENDULUM. 




(1) 




4 



limit of small driving amplitude a/I < 1 and describe how the dynamical behavior qualitatively 
changes as a function of the driving frequency. 

For a small amplitude drive the lower equilibrium at 9 = remains stable unless particular 
parametric resonance conditions ^ ~ ^ with n = integer, are me1>^. As we increase the frequency 
of the external drive from 7 ~ 2loq we observe qualitatively different regimes. First the motion in 
phase space is completely chaotic and both the lower and upper equilibrium {9 = 0, ir) are unstable, 
then the lower equilibrium becomes stable while the upper equilibrium remains unstable, finally 
when f J- > y/2 both the upper and lower equilibrium are stable. The surprising phenomenon that 
the upper position becomes stable (and the pendulum performs oscillations around this inverted 
position) is known in the literature as dynamical stabilization and was first explained by Kapitza. 
He showed that for small amplitude and high frequency driving the dynamic of the driven pendulum 
can be accurately described by a time-independent effective Hamiltonian, moreover the effective 
potential energy develops a local minimum at 9 = ir when f ^ > V% explaining the oscillations 
around the inverted position. 

Usually the dynamical stabilization is obtained by splitting the degrees of freedom into fast and 
slow modes, eliminating the fast modes, and obtaining the effective potential for the slow modest. 
This procedure has a limitation that it can not be easily extended to either interacting systems 
or to the quantum domain. It is also unclear whether averaging over fast degrees of freedom will 
lead to the Hamiltonian equation of motion in each order of the expansion. Here we show that 
the dynamical stabilization phenomenon can be understood through the Magnus expansion of the 
quantum evolution operator in powers of the inverse frequency. The advantage of this method is 
that allows us to analyze behavior of the periodically driven interacting systems. 

Magnus Expansion. Let us now start by reviewing the basis aspects of the Magnus Expan- 



sion^ (ME) (for an excellent review see Ref. 20]). This expansion equally applies to quantum and 
classical systems, but it is more convenient to use the language of quantum mechanics. We imagine 
that the system is prepared initially in some quantum-mechanical state |Y>o), which evolves in time 
under some periodic Hamiltonian H(t) with a period T. Then if we look stroboscopically into the 
system at times t n = nT its wave function will be given by 

|^(nT)) = U{TT\^) (2) 

where U(T) is the evolution operator during one period of oscillation: 

U(T) = T T exp (~ £ dTH{t)^j . (3) 
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Here the time ordering exponent implies that the later times in the integral always appear on the 
left. Since the evolution operator is unitary it can be always represented as the exponent of some 
Hermitian operator, which we can identify with the effective Floquet Hamiltonian H e ff. 

U(T) = exp (~H eff T) (4) 

The problem of finding the Floquet Hamiltonian is thus equivalent to the problem of evaluating 
the logarithm of a time-ordered exponential. For some simple noninteracting Hamiltonians it is 
actually possible to find the Floquet Hamiltonian explicitly and this was recently used to predict 
new topological states of light and atom o 21 ' 22 . But in general it is a very hard problem with no 
closed form solution because the Hamiltonian at different times does not commute with itself. 
Unlike in statistical physics this problem remains equally hard even in the classical limit because 
of the singular factor 1/h in the exponent. 

A possible root to compute the Floquet Hamiltonian is to expand the evolution operator in 
Taylor series in the period T and then take the logarithm of this Taylor series. The corresponding 
M E 19 ' 20 is guaranteed to converge for finite-dimensional Hamiltonians and sufficiently short periods. 
We briefly review this expansion and explicitly compute few terms for a specific models in the 
Supplementary Material. In thermodynamic limit the convergence of the ME can not be proven in 
general. This situation is very similar to the high-temperature expansion in statistical physics 2 ^. 
Usually the breakdown of the high temperature expansion is an indication of a phase transition. 
For simple few body systems there are many physical situations where the ME converges quickly 
to the exact Floquet Hamiltonian and only few terms are necessary to describe the dynamics. For 
this reason the ME found many applications in different areas of physics and mathematics (see the 
review and refs. therein). The first two terms in the ME read: 

H eff = T So dt H(h), ilfjf = 2rpy n o< t 2< t 2 <T dt i dt 2 H(ti),H(t 2 ) • (5) 
Higher order terms have a similar structure containing higher order commutators. They are mul- 
tiplied by higher powers of the period and hence the period of the driving plays the role similar to 
the inverse temperature in the high temperature expansion. For classical systems the equivalent 
expansion can be obtained by substituting the commutators between the operators with the Pois- 
son brackets: [. . . ]/ih — > {•••}• In the case of the Kapitza pendulum it is sufficient to compute the 

eff + a eff + n eff 



first three terms in the ME (see Supplemental Material), i.e. Heff ~ Hjiff + ^iff ^eff They 



read: 

h eff~2^^ -" la, LU ^. "eff 



H^ f = ^P 2 -m^ cose, H$ f = 



( c \ 

H?f)=r^+™i({^Y-i)^ 2 e 



We point that up to this order, the ME is the same for quantum and classical pendulum. Here, 
the first term is simply the time averaged Hamiltonian, which reduces to the Hamiltonian of the 
undriven pendulum. The second term is zero. Actually it can be shown that for symmetric pulses 
all even terms in the ME are identically equal to zer o 24 ' 25 . The third term H^X is the first nontrivial 
contribution to the Floquet Hamiltonian, which is sufficient to explain the dynamical stabilization 
of the pendulum found by Kapitza. Indeed in the limit a <C I the Floquet Hamiltonian describes 
a pendulum moving in the effective potential U e ff = mujQ ^— cos 9 + ^2^) sm2 ■ ^ ^ s eas y to 
see that this potential develops a new minimum at the inverted position, 8 = ir when the condition 
j ^ > y/2 is met. Moreover one can check that the higher order terms in the ME remain small as 
long as this condition is met. When the period of the external driving is increased (the frequency 
is reduced) higher order terms in the ME need to be computed to accurately describe the dynamics 
of the system. If the period is increased even further there is no guarantee that the ME converges 
to the effective Floquet Hamiltonian. In this regime we observe that the motion of the pendulum 
becomes chaotio^. 



III. ENERGY LOCALIZATION TRANSITION IN INTERACTING SPIN CHAINS. 

In the rest of the paper we will focus an Id interacting classical and quantum spin models 
with periodic boundary conditions which is driven by periodically quenching between two different 
Hamiltonians, Hq and H\. Although we believe that our results are general, we prefer to keep the 
discussion focused on this specific system. In Fig. [2] we show two equivalent descriptions of the 
same protocol. The left panel represents the actual time dependent sequence of pulses switching 
between the two Hamiltonians. On the right panel we depict an equivalent protocol where instead 
we are performing a single quench to the Floquet Hamiltonian and then a quench back to the 
original Hamiltonian Hq at the time of measurement. The second quench is not necessary, it 
highlights the fact that we do the measurements of observables like energy with respect to the 
Hamiltonian Hq. 

Without loss of generality we can choose the Hamiltonian Hq to be a simple precession in the 
external magnetic field and the Hamiltonian H\ to be interacting and ergodic: 
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FIG. 2: Two equivalent description of the driving protocol: (left) sequence of sudden quenches between Ho 
and H\ and (right) single quench from Hq to the effective Floquet Hamiltonian H e ff and back to Hq. 



where, we have defined the shorthand notations:: 



Hbh 



Yn S ni Hz — Yn i S n s n+l) ) H\\ ~ Yn { S n s n+1 + s ™ S n+l) 
H'z = Yn \ s n s n+2) ' H\\ = Yi { s n s n.+2 + s " s n+2) 

Let us point that this system is invariant under space translation and ir — rotation around the 
x — axis (s* — > s^, s„ — > —s n , — > —s^). For numerical calculations we choose the following 
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j. We checked that our results are not tied 



parameters: B x = 1, J z = — Jl 
to any particular choice of couplings. 

As pointed out earlier, we can expect two qualitatively different regimes depending on the 
period of the driving. At long periods the system has enough time to relax to the stationary state 
between the pulses and thus is expected to constantly absorb energy until it reaches the infinite 
temperature. This situation is similar to what happens for driving with random periods^. On the 
contrary if the period is very short we can expect that the Floquet Hamiltonian converges to the 
time averaged Hamiltonian. Since the whole time evolution can be viewed as a single quench to 
the Floquet Hamiltonian (right panel in Fig. [2]) we expect that the energy will be localized even in 
the infinite time limit as long as the Floquet Hamiltonian is well defined and local. Noticing that 
the commutator of two local extensive operators is local and extensive we see from Eq. ([5]) that 
the Floquet Hamiltonian is local an extensive in each order of ME. Thus the question of whether 
the energy of the system is localized in the infinite time limit or reaches the maximum possible 
value is tied to the question of convergence of the ME. To our knowledge there are no statements 
in the literature about this convergence in the thermodynamic limit. Therefore we will rely on the 
specific spin model to establish that for short driving periods the ME indeed converges and the 
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energy of the system is localized while for longer periods this expansion diverges and the systems 
is heated towards infinite temperature. 

For the quantum system we consider a spin-| chain initially polarized along the x-axis, which is 
the ground state of Hq. We then compute the time evolution operator by exact diagonalization of 
the two Hamiltonians Hq and H\. Translational invariance of the system allows us to restrict the 
analysis to the zero momentum sector and analyze chains with the number of spins up to N = 17, 
which contains 7712 states. For the system of classical spins we numerically integrate the Bloch 
equations derived from the Hamiltonians (J7J) using the forth-order Runge-Kutta algorithm. We took 
the initial conditions which maximally mimic the Wigner function of the quantum chain, namely 
s x = 2 ' s v = cos 01 anc ^ Sz = \72 s * n a w ^ ere < a < 2-7T is a uniformly distributed random 
variable. The advantage of quantum systems is that it is straightforward to exactly obtain the 
infinite time limit of the evolution by projecting the initial state to the eigenstates of the effective 
Hamiltonian, while the advantage of the classical systems is that one can analyze significantly 
larger system sizes. 

In Fig. [3] we show the excess energy in the infinite time as a function of the pulse times To and 
T\. The blue (orange) regions in parameter space correspond to the system being unable (able) 




TO 

FIG. 3: (Color online) Excess energy of the quantum spin chain in the long time limit: Q = 
(ip(t)\Ho\tp(t)}t^oo ~ E gs , where E gs is the ground state energy of the Hamiltonian Hq, as a function 
of the pulse times To and T\ in units of h/B x . Blue (orange) regions correspond to small (large) excess 
energy. The data is obtained by the exact diagonalization of a spin-i chain with N = 15 spins. 

to absorb energy from the time-dependent driving. The periodicity in To, which is clearly visible 
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in Fig. [3l is due to the periodicity of the spin precession generated by the Hamiltonian Hq. For 
this reason we focus only on To < ix. This plot illustrates a sharp crossover between localized 
and delocalized phases of the driven spin chain as a function of the pulse times. To establish that 
this crossover becomes a true phase transition in the thermodynamic limit we next show detailed 
analysis of the excess energy along two line-cuts in the T$ — T\ plane. Close to the To — axis (green 
arrow in Fig. [3]), the effective Hamiltonian H e ff can be computed to the first order in T\ and all 
orders in To using the Hausdorff-Baker-Campbell formula 2 ^. For the present problem it is possible 
to obtain the result of this resummation in the closed analytic form (see Methods): 



where H av = j, (HqTq + H\Ti) is the time-averaged Hamiltonian, M is a local operator that couples 
nearest-neighbor and next-nearest-neighbor spins and A = Bx ^° ■ 

From Eq. ([8]) we see that the effective Hamiltonian becomes singular for To = nil (in units 
of h/B x) no matter how small T\ or the coupling constants in H\ (the Js) are. The location of 
the singularity is also manifestly independent of the system size. As it is seen from the numerical 
simulations (Fig. [3]), the singularity in the effective Hamiltonian is also manifested in the excess 
energy of the system in the infinite time limit. We point out that, in the limit of small Tl, our 
system directly extends the kicked rotor model to the many-spin domain. Indeed most of the time 
the spins precess around the magnetic field B x getting periodically short kicks by the interacting 
Hamiltonian H\. Thus in this limit the many-body localization transition directly generalizes the 
well known kicked rotor localization transitio n 11 ' 12 . 

Away from the To — axis, nested commutators of order T" for n > 1 need to be included in 
the effective Hamiltonian. These commutators become difficult to compute analytically at high n. 
Like in the high temperature expansion in statistical physics they involve multiple spin interactions 
and become non-local in space. Therefore we have to rely on numerics. In Fig. [3] we analyze the 
long time limit of the excess energy along the generic direction T\ = Tq — 2 for 2 < To < 3 (pink 
arrow in Fig. [3]). 

In Fig. [5] we show the error between the exact asymptotic value of the normalized excess energy 
(see Fig. H|) and the corresponding excess energies obtained by truncating the ME. From the left 
panel we see that, as expected, for short periods the ME becomes asymptotically exact. From the 
right panel we see that the truncation error at short periods does not increase with the system 
size indicating existence of the localized phase in the thermodynamic limit. The latter implies 
existence of the localization transition. For long periods (see Fig. [4]) the energy clearly approaches 
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FIG. 4: (Color online) Asymptotic value of the normalized excess energy, 2 ^ H °^ t ~°^ E E - ) 9S where E gs and 
Emax are the lowest and highest eigenvalues of Hq. The exact results for different system sizes arc compared 
with the predictions obtained by truncating the Magnus Expansion to different orders: H e ff(k) contains 
terms of order Tq™T™ with m + n < k and -ff e //(oo + 5) denotes the non-perturbative result Eq. ((5J together 
with all the terms Tq 1 T™ with m + n < 5. The results from the Magnus Expansion for different system sizes 
are identical within the image resolution. 

the infinite temperature asymptotic value. 

The localization transition is also found for the system of classical spins, see Fig. EJ On the 
left panel we show the time evolution of the excess energy for different values of the pulse times 
To and T\. We observe that close to the transition the evolution of the excess energy is extremely 
slow (note the scale in Fig. ED suggesting that the this transition is characterized by a divergent 
time scale, i.e. the time required to reach the asymptotic state. On the right panel we show the 
asymptotic value of the excess energy for different system sizes. These asymptotic values have been 
computed by averaging the energy over the last 10 cycles of the evolution. However close to the 
transition, it is not clear if the energy has reached its asymptotic value. So more work is needed 
to check whether the sharp crossover in classical systems seen on the right panel of Fig. [6] becomes 
a true phase transition in the long time limit. The analytical argument based on the resummation 
of the Hausdorff-Baker-Campbell formula (Eq. flSJ)) suggests this is in fact a transition. 
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FIG. 5: (Color online) Difference between the asymptotic value of the exact normalized excess energy and 
the approximate values obtained from the truncated Floquet Hamiltonians at different orders, for a fixed 
system size N = 15 (left panel). Same error for fixed order of truncation versus system size (right panel). 
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FIG. 6: (Color online) Time evolution of the excess energy for spin, Q = > f° r a classical chain with 

N = 128 spins (left panel) and asymptotic excess energy for different systems sizes (right panel). 



IV. CONCLUSIONS. 



We presented analytic and numerical arguments suggesting existence of the localization transi- 
tion in interacting driven systems in the thermodynamic limit. For short periods of driving T < T c 
the system is unable to absorb energy beyond some threshold even in the infinite time limit while 
for longer periods T > T c the system absorbs energy and becomes delocalized in the entire Hilbert 
space. This interpretation is confirmed by the study of the diagonal entropy^ (see Supp. Materi- 
als): Sd{t) = — J2kPo(t) logPoW wnere Po(*) * s * ne occupation probability at time t of the k — th 
eigenstate of Hq. The diagonal entropy serves as a measure of the occupation of the Hilbert space. 
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For short periods it remains bounded at all times and at long periods, up to small sub-extensive 
corrections, approaches the infinite temperature entropy given by the logarithm of the total number 
of states. We associate the derealization transition with the divergence of the short time ME for 
the effective Floquet Hamiltonian similarly to the divergence of the high-temperature expansion in 
statistical physics. In this sense this transition is reminiscent of the divergence of the short time 
expansion in the Loschmidt echo recently obtained for a quench in the transverse field Ising spin 
chain22. More work is needed to precisely characterize the two phases, the nature of the transition 
and the precise connection of the transition with the radius of convergence of the MFj2£. While at 
the moment we can not rigorously prove that this is a true phase transition, the right panel in Fig. [5] 
gives a very convincing evidence that at short periods the energy always remains bounded even 
in the thermodynamic limit indicating the existence of the localized phase. We believe that this 
transition is generic for systems with bounded single particle excitations. This transition should be 
experimentally detectable in isolated interacting systems like cold atoms, cold ions, nuclear spins 
and various materials with interacting spin degrees of freedom. 
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V. METHODS. 

Here, we show how to resum the Hausdorff-Baker-Campbell formula (HBC)— to obtain the 
Hamiltonian ([8]). First, it is convenient to introduce the shorthand notations X — > —-^HqTq, Y — > 
— jrHiTi and Z — > — jK e jjT where T = Tq + T\. With these notations the effective Hamiltonian 
for the protocol considered (see Fig. [2]) is defined by Z = log [e x e Y e x ~\ . Expanding this expression 
to the first order in Y we obtain the HBC formula: 

°° o/o2n— 1 _ i \ 

Z = 2X + Y - (2n) , B 2n {X 2n ; Y) + 0(Y 2 ) (9) 

n=l ^ ' 

where are the first Bernoulli numbers (Bq = 1,B\ = —^,B^ = B§ = ... = 0) and 
Y : [X; ... [X; Y]]] is the nested commutator with k operators X and only one operator Y appearing 



13 



at the right most position. Note that Eq. Q is different from the expression usually reported in 
the literature, where the less symmetric situation Z = log [e x e Y ~\ is considered. The simple form 
of the Hamiltonian Hq allows us to compute all the nested commutator analytically. In fact, using 
the explicit form of Hq and H\ (see Eq. ([7])) together with the commutation relations for spin 
operators [s"; Sm] = i8 njm e a pjS% , it is easy to verify that 

{{H Bx ) 2n ;H z } = 2^W zz - yy , {(H Bx ) 2n ;H' z ) = 2^W' zz _ yy 
{(H Bx ) 2n ;H {l ) = -2^W zz . yy , {(H Bx ) 2n ;H' {l } = -2^W' zz _ yy 

where W zz - yy = E„ (<-i< " s Li s ») and w L-yy = E„ (<-2< ~ S n-2 S «) • Plugging these 
expressions in Eq. and performing the sum we find Eq. ([8]) where M = (<7 Z — Jm) W zz - yy + 
(^J' z — J|'^ W zz _ yy . We note that for spin rotation invariant situations, J z = Jy and J' z = Jjj, M = 
and the effective Hamiltonian reduces to the time-averaged Hamiltonian. This is expected since in 
this case Hq and Hi commute with each other. 
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I. MAGNUS EXPANSION FOR THE KAPITZA PENDULUM 



Here we will show how to apply the Magnus Expansion (ME) to the Kapitza pendulum problem. 
In particular we will show how to derive the dynamical stabilization condition, first obtained by 
Kapitza. The quantum time-dependent Hamiltonian for the Kapitza pendulum is: 

1 



H(t) = -^pj + f(t) cos 
2m 



(1) 



where f(t) 
relations 
article^) : 



-m [uj, 



We 



o + f 7 2 cos (7*)) an d Pe are quantum operators with canonical commutation 
ih. The explicit form of the first three terms in the ME are (see the review 
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where the time integration domains are ordered, i.e. < t n < t n -\ < ... < t\ <T. Recalling that 



the period of the driving is T = ^ after some simple algebra we obtain: 
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Substituting the explicit value for the commutators in H pf f we obtain: 
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Combining Eqs. (J3j> an d Eq. (JH) we obtain the first three terms in ME for the quantum Kapitza 
pendulum. Up to this order, the classical ME (see Eq. (6) in the main text) can be obtained from 
the quantum counterpart by substituting the quantum operators with classical variables. This is 
not true in general and a more rigorous approach is necessary to derive the classical limit of the 
ME. 

Before showing the general approach to obtain the classical limit of the ME we note that the 
first three terms in the ME suffice to explain the dynamical stabilization of the classical Kapitza 

2 2 2 

pendulum. Let us assume that -j- <C -^r- (we will check this assumption a posteriori) then by 

rW and 



collecting the terms in H^X and H^fX that involves only the coordinates we obtain the effective 



potential U e ff = ttiuq — cos 9+ ^2^) sin 2 which develops a new minimum at the inverted 
position, 6 = ir, for f ^- > \/2. Moreover close to the stabilization transition {^j-) 2 ~ 2a; 2 and we 
can check that our assumption is justified as long as j <C ^, which is one of the condition required 
by the original derivation by Kapitza. 

Now let us discuss how one can obtain obtain the classical limit of the ME through the phase- 
representation of Quantum Mechanics. Let us briefly review the formalism^. In the phase-space 
representation of quantum mechanics the quantum operator Ci is replaced by its Weyl Symbol 
which is a classical function over the phase-space variables x and p (which can be vectors for 
multi-dimensional problems). The Weyl symbol is defined as 



/s * s 
ds (x - -\&\x + -) exp 



h P ' S 



(5) 



If the quantum operator £l(x,p) is written in the symmetrized form then it can be shown that its 
Weyl symbol £l w is simply obtained by the substitution x — > x and p — > p. In particular, this is 
true for all operators of the form £l(x,p) = A{x) + B(p). In this representation the commutator of 
two quantum operators is written in terms of the Moyal Bracket: 

Cli;Cl 2 ->• ih{i}i, w ;£l 2 ,w} MB = ih £li )W ^-^sin f^O^ Q 2 ,w (6) 

where A is the differential operators A = c^J d x — d p . By expanding —2 sin (f A) in powers of h 
we obtain the classical limit, i.e. the Poisson Bracket (zero order in h), and the quantum correction 
(higher powers of h). Note that the first quantum correction is proportional to A 3 . Thus if we 
are interested in the classical limit of ME {h — > 0) it suffices to truncate the expansion to order 
A. This amount to the replacement of the quantum commutator with the Poisson Bracket, i.e. 
4J...] — > {...}. In the case of Kapitza pendulum, the only non-zero contribution to the first few 
nested commutators comes from the order A and the classical and quantum ME are identical. 
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This is not always true since at higher order in the ME terms of order A 3 (or higher) need to 
be considered and non-trivial differences between the classical and quantum expansions will be 
generated. 



II. EXPLICIT FORM OF FIRST TERMS IN THE MAGNUS EXPANSION FOR THE 

INTERACTING QUANTUM SPIN CHAIN 

As mentioned in the main text, all the even order terms in the ME vanish for a time-symmetric 
driving 3 -^, i.e. H(t) = H(T — t). Here, for completeness, we report the first not-vanishing terms 
in the ME for the protocol considered in the main paper. First, it is convenient to introduce the 
shorthand notations X — > —-^HqTq, Y — > — iHiTi and Z — > —jrH e ffT where T = Tq + T\. Then 
by the definition of evolution operator, the effective Hamiltonian is defined by Z = log [f7(T)] = 



log [e x e Y e x ~\ . Using a Mathematica code similar to the one described in Ref. 
computed the first few terms in the ME: 

Z = (2X + Y) + \ (- {XXY} - {YXY}) + 
b 



we have explicitly 



-t-[7{XXXXY} + {YYYXY} + 6{XYXXY} (7) 
+ 8{YXXXY} + 12 {YY XXY} - 4{XYYXY}) + ... 



where the symbol {XYXXY} is a short hand notation for the right nested commutators 
[X, [Y, [X, [X,Y]]]]. We note that there is no unique way of writing Z in terms of nested com- 
mutators due to the Jacoby identity of nested commutators. For example, it is easy to show 
that {YYXXY} = {YXYXY}. Substituting the definition for X,Y,Z in Eq. © we obtain 
Heff = II, j f + ir] f + "If + - where: 

H% = i (H T + HiTx) 

A3) _ 1 ((T \Z 



H eff = WT ((f ) T i {HoHoH,} + (f ) If {H.HvH,}^ 
H eff = 360W ({?£)' l T 1 7{H HoH HoH 1 }+(^)T?{H 1 H 1 H 1 H H 1 }+ (8) 
(f ) 3 If (6 {HoHxHoHqHx} + 8 {#i# #o#o#i}) + 
(^) 2 T 3 (12 {HxHxHoHoHx} - 4 {HqH^HqHx}) 



These terms are different from the ones usually reported in the literature where the less symmetric 
protocol, Z = log [e x e Y ~\ , is considered. 
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III. BEHAVIOR OF THE DIAGONAL ENTROPY ACROSS THE LOCALIZATION 
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FIG. 1: (Color online) (Color online) Asymptotic value of the normalized diagonal entropy, Sd/S max where 
Smax is the logarithm of the total number of the microstates. The exact results for different system sizes are 
compared with the predictions obtained by truncating the Magnus Expansion to different orders: H e ff{k) 
contains terms of order Tq™T" with m + n < k and iJ e //(oo + 5) denotes the non-pcrturbative result in To 
(Eq. (8) in the main text) together with all the other terms T™T™ with m + n < 5. The results from the 
Magnus Expansion for different system sizes arc identical within the image resolution. 

In this Section we show our numerical results for the behavior of the diagonal entropy across the 
localization transition. The diagonal entropy Q] is defined as Sd(t) = — Po(t) logp^t) where 
Pg'(t) is the occupation probability at time t of the k — th eigenstate of Hq. The diagonal entropy 
serves as a measure of the occupation in the Hilbert space. As in Fig. 4 in the main text we study 
the behavior along the generic line T\ = Tq — 2 for 2 < To < 3. In particular, in Fig. Q] we show the 
asymptotic value the diagonal entropy for different system sizes and we compare it to the values 
obtained by truncating the ME to different orders. The asymptotic values have been computed 
by projecting the initial state to the eigenstates of H e ff and then back to Hq. This procedure is 
equivalent to the assumption of infinite time averaging with respect to H e ff and the asymptotic 
values obtained correspond to the prediction of the diagonal ensemble of H e ff. Fig. [T] shows that 
the energy increase observed in Fig. 4 of the main text is indeed caused by a derealization of 
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the system in Hilbert space. Finally we note that for short (long) periods the normalized entropy 
decreases (increases) with the system size. This suggests that the diagonal entropy might be a 
better indicator than the excess energy to precisely locate the transition point. 



IV. EVOLUTION OF ENERGY AND ENTROPY FOR THE INTERACTING 

QUANTUM SPIN CHAIN 
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FIG. 2: (Color online) Time evolution of (left) the excess energy for spin Q = 2 — ^ 

Emax are the lowest and highest eigenvalues of Hq and (right) the normalized diagonal entropy, Sd(t) / Smax 
where Smax is the logarithm of the number of states in the system. The dotted lines indicate the predictions 
from the diagonal ensemble. 



In this section we report our numerical findings for the time-evolution of the excess energy and 
the diagonal entropy of the quantum chain with N = 17 spins for different values of the pulse times 
Tq = T\ = T/2. In particular, in Fig. [2] we compare the time series of these quantities with the 
predictions from the diagonal ensemble of H e f f . 

On the left panel we observe that the excess energy, Q = 2 ^• H °^ t ~°^ E E ^ a , approaches the 
diagonal ensemble predictions for all the values of the pulsed times To and T\ . Moreover, we see that 
for short and long periods (such as Tq = T\ = 1.02 and Tq = T\ = 2.3 respectively) the asymptotic 
value of the excess energy is approached quickly while for intermediate periods (such as Tq = T\ = 
1.92) this approach is extremely slow. This seems to suggest that the localization transition is 
associated with a divergent time scale, i.e. the time to approach the asymptotic value. On the 
right panel we observe a similar behavior for the diagonal entropy, Sd(t) = — YlkPoft) l°SPo(*) 
where Pq(^) ^ s the occupation probability at time t of the k — th eigenstate of Hq. In this case 



however the prediction from the diagonal ensemble seems to over-estimate, for long periods, the 
value of the diagonal entropy. This small discrepancy, which is non-extensive and thus vanishes 
in thermodynamic limit^, is caused by the fact that the entropy is not a linear function of the 
density matrix. Finally we note that for long periods the fluctuations in the time-series of the both 
observables (excess energy and diagonal entropy) are suppressed. This is expected since in the 
delocalized phase more states participate in the dynamics and the statistical average is improved. 
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